Determinants of immunoglobulin G responses to respiratory syncytial virus and rhinovirus in children and adults.

Introduction: Exposure to respiratory viruses is a significant cause of morbidity and affects virus-specific antibody levels. Little is known about determinants associated with immune response to these viruses. We aimed to investigate the determinants of respiratory syncytial virus (RSV)- and rhinovirus (RV)- specific IgG responses in both children and adults. Methods: The study is based on the EGEA cohort, composed of 530 samples of children in EGEA1 (1991-95) and 1241 samples of adults in EGEA2 (2003-07). Cumulative RV-specific IgG levels (species A, B and C) and IgG levels to RSV-G protein were measured by using micro-array technoloy. Multiple linear mixed models (random effect to account for familial dependence) were performed to assess associations between age, sex, body mass index (BMI), tobacco smoke exposure and season of blood sampling with RSV-and RV-specific IgG levels. Results: In children (11.1 ± 2.8 years old, 57% boys), higher RV-specific IgG levels were associated with older age (only for RV-B), female sex and lower BMI, while only older age was associated with higher RSV-specific IgG levels. In adults (43.5 ± 16.7 years old, 48% men), younger age, female sex, lower BMI, active smoking and all seasons except summer were associated with higher RV-specific IgG levels. Older age, active smoking and all seasons except summer were associated with higher RSV-specific IgG levels. Conclusion: Personal and seasonal determinants of RSV- and RV-specific IgG levels seem to vary according to the respiratory virus type and between children and adults, suggesting different patterns of responses along the life course.

Network propagation for GWAS analysis: a practical guide to leveraging molecular networks for disease gene discovery.

MOTIVATION: Genome-wide association studies (GWAS) have enabled large-scale analysis of the role of genetic variants in human disease. Despite impressive methodological advances, subsequent clinical interpretation and application remains challenging when GWAS suffer from a lack of statistical power. In recent years, however, the use of information diffusion algorithms with molecular networks has led to fruitful insights on disease genes. RESULTS: We present an overview of the design choices and pitfalls that prove crucial in the application of network propagation methods to GWAS summary statistics. We highlight general trends from the literature, and present benchmark experiments to expand on these insights selecting as case study three diseases and five molecular networks. We verify that the use of gene-level scores based on GWAS P-values offers advantages over the selection of a set of 'seed' disease genes not weighted by the associated P-values if the GWAS summary statistics are of sufficient quality. Beyond that, the size and the density of the networks prove to be important factors for consideration. Finally, we explore several ensemble methods and show that combining multiple networks may improve the network propagation approach.

TNS1 and NRXN1 Genes Interacting With Early-Life Smoking Exposure in Asthma-Plus-Eczema Susceptibility.

PURPOSE: Numerous genes have been associated with allergic diseases (asthma, allergic rhinitis, and eczema), but they explain only part of their heritability. This is partly because most previous studies ignored complex mechanisms such as gene-environment (G-E) interactions and complex phenotypes such as co-morbidity. However, it was recently evidenced that the co-morbidity of asthma-plus-eczema appears as a sub-entity depending on specific genetic factors. Besides, evidence also suggest that gene-by-early life environmental tobacco smoke (ETS) exposure interactions play a role in asthma, but were never investigated for asthma-plus-eczema. To identify genetic variants interacting with ETS exposure that influence asthma-plus-eczema susceptibility. METHODS: To conduct a genome-wide interaction study (GWIS) of asthma-plus-eczema according to ETS exposure, we applied a 2-stage strategy with a first selection of single nucleotide polymorphisms (SNPs) from genome-wide association meta-analysis to be tested at a second stage by interaction meta-analysis. All meta-analyses were conducted across 4 studies including a total of 5,516 European-ancestry individuals, of whom 1,164 had both asthma and eczema. RESULTS: Two SNPs showed significant interactions with ETS exposure. They were located in 2 genes, NRXN1 (2p16) and TNS1 (2q35), never reported associated and/or interacting with ETS exposure for asthma, eczema or more generally for allergic diseases. TNS1 is a promising candidate gene because of its link to lung and skin diseases with possible interactive effect with tobacco smoke exposure. CONCLUSIONS: This first GWIS of asthma-plus-eczema with ETS exposure underlines the importance of studying sub-phenotypes such as co-morbidities as well as G-E interactions to detect new susceptibility genes.

Trait selection strategy in multi-trait GWAS: Boosting SNPs discoverability.

Since the first Genome-Wide Association Studies (GWAS), thousands of variant-trait associations have been discovered. However, the sample size required to detect additional variants using standard univariate association screening is increasingly prohibitive. Multi-trait GWAS offers a relevant alternative: it can improve statistical power and lead to new insights about gene function and the joint genetic architecture of human phenotypes. Although many methodological hurdles of multi-trait testing have been discussed, the strategy to select trait, among overwhelming possibilities, has been overlooked. In this study, we conducted extensive multi-trait tests using JASS (Joint Analysis of Summary Statistics) and assessed which genetic features of the analysed sets were associated with an increased detection of variants as compared to univariate screening. Our analyses identified multiple factors associated with the gain in the association detection in multi-trait tests. Together, these factors of the analysed sets are predictive of the gain of the multi-trait test (Pearson's ρ equal to 0.43 between the observed and predicted gain, P < 1.6 × 10-60). Applying an alternative multi-trait approach (MTAG, multi-trait analysis of GWAS), we found that in most scenarios but particularly those with larger numbers of traits, JASS outperformed MTAG. Finally, we benchmark several strategies to select set of traits including the prevalent strategy of selecting clinically similar traits, which systematically underperformed selecting clinically heterogenous traits or selecting sets that issued from our data-driven models. This work provides a unique picture of the determinant of multi-trait GWAS statistical power and outline practical strategies for multi-trait testing.

Identification of novel genes influencing eosinophil-specific protein levels in asthma families.

BACKGROUND: Eosinophils play a key role in the asthma allergic response by releasing cytotoxic molecules such as eosinophil cationic protein (ECP) and eosinophil-derived neurotoxin (EDN) that generate epithelium damages. OBJECTIVE: We sought to identify genetic variants influencing ECP and EDN levels in asthma-ascertained families. METHODS: We performed univariate and bivariate genome-wide association analyses of ECP and EDN levels in 1018 subjects from the EGEA study with follow-up in 153 subjects from the Saguenay-Lac-Saint-Jean study and combined the results of these 2 studies through meta-analysis. We then conducted Bayesian statistical fine mapping together with quantitative trait locus and functional annotation analyses to identify the most likely functional genetic variants and candidate genes. RESULTS: We identified 5 genome-wide significant loci (P &lt; 5 × 10^-8) including 7 distinct signals associated with ECP and/or EDN levels. The genes targeted by our fine mapping and functional search include RNASE2 and RNASE3 (14q11), which encode EDN and ECP, respectively, and 4 other genes that regulate ECP and EDN levels. These 4 genes were JAK1 (1p31), a transcription factor that plays a key role in the immune response and acts as a potential therapeutic target for eosinophilic asthma; ARHGAP25 (2p13), which is involved in leukocyte recruitment to inflammatory sites; NDUFA4 (7p21), which encodes a component of the mitochondrial respiratory chain and is involved in cellular response to stress; and CTSL (9q22), which is involved in immune response, extracellular remodeling, and allergic inflammation. CONCLUSION: Analysis of specific phenotypes produced by eosinophils allows the identification of genes that play a major role in allergic response and inflammation, and offers potential therapeutic targets for asthma.

Genome-Wide Association Study of Fluorescent Oxidation Products Accounting for Tobacco Smoking Status in Adults from the French EGEA Study.

Oxidative stress (OS) is the main pathophysiological mechanism involved in several chronic diseases, including asthma. Fluorescent oxidation products (FlOPs), a global biomarker of damage due to OS, is of growing interest in epidemiological studies. We conducted a genome-wide association study (GWAS) of the FlOPs level in 1216 adults from the case-control and family-based EGEA study (mean age 43 years old, 51% women, and 23% current smokers) to identify genetic variants associated with FlOPs. The GWAS was first conducted in the whole sample and then stratified according to smoking status, the main exogenous source of reactive oxygen species. Among the top genetic variants identified by the three GWAS, those located in BMP6 (p = 3 × 10-6), near BMPER (p = 9 × 10-6), in GABRG3 (p = 4 × 10-7), and near ATG5 (p = 2 × 10-9) are the most relevant because of both their link to biological pathways related to OS and their association with several chronic diseases for which the role of OS in their pathophysiology has been pointed out. BMP6 and BMPER are of particular interest due to their involvement in the same biological pathways related to OS and their functional interaction. To conclude, this study, which is the first GWAS of FlOPs, provides new insights into the pathophysiology of chronic OS-related diseases.

Analysis of GWAS-nominated loci for lung cancer and COPD revealed a new asthma locus.

BACKGROUND: Asthma, lung cancer (LC) and chronic obstructive pulmonary disease (COPD) are three respiratory diseases characterized by complex mechanisms underlying and genetic predispositions, with asthma having the highest calculated heritability. Despite efforts deployed in the last decades, only a small part of its heritability has been elucidated. It was hypothesized that shared genetic factors by these three diseases could help identify new asthma loci. METHODS: GWAS-nominated LC and COPD loci were selected among studies performed in Caucasian cohorts using the GWAS Catalog. Genetic analyses were carried out for these loci in the Saguenay-Lac-Saint-Jean (SLSJ) asthma familial cohort and then replicated in two independent cohorts (the Canadian Cohort Obstructive Lung Disease [CanCOLD] and the Epidemiological Study of the Genetics and Environment of Asthma [EGEA]). RESULTS: Analyses in the SLSJ cohort identified 2851 and 4702 genetic variants to be replicated in the CanCOLD and EGEA cohorts for LC and COPD loci respectively. Replication and meta-analyses allowed the association of one new locus with asthma, 2p24.3, from COPD studies. None was associated from LC studies reported. CONCLUSIONS: The approach used in this study contributed to better understand the heritability of asthma with shared genetic backgrounds of respiratory diseases.

Identification of OCA2 as a novel locus for the co-morbidity of asthma-plus-eczema.

BACKGROUND: Numerous genes have been associated with the three most common allergic diseases (asthma, allergic rhinitis or eczema) but these genes explain only a part of the heritability. In the vast majority of genetic studies, complex phenotypes such as co-morbidity of two of these diseases, have not been considered. This may partly explain missing heritability. OBJECTIVE: To identify genetic variants specifically associated with the co-morbidity of asthma-plus-eczema. METHODS: We first conducted a meta-analysis of four GWAS (Genome-Wide Association Study) of the combined asthma-plus-eczema phenotype (total of 8807 European-ancestry subjects of whom 1208 subjects had both asthma and eczema). To assess whether the association with SNP(s) was specific to the co-morbidity, we also conducted a meta-analysis of homogeneity test of association according to disease status ("asthma-plus-eczema" vs. the presence of only one disease "asthma only or eczema only"). We then used a joint test by combining the two test statistics from the co-morbidity-SNP association and the phenotypic heterogeneity of SNP effect meta-analyses. RESULTS: Seven SNPs were detected for specific association to the asthma-plus-eczema co-morbidity, two with significant and five with suggestive evidence using the joint test after correction for multiple testing. The two significant SNPs are located in the OCA2 gene (Oculocutaneous Albinism II), a new locus never detected for significant evidence of association with any allergic disease. This gene is a promising candidate gene, because of its link to skin and lung diseases, and to epithelial barrier and immune mechanisms. CONCLUSION: Our study underlines the importance of studying sub-phenotypes as co-morbidities to detect new susceptibility genes.

[Tobacco and DNA methylation: the case for epigenetic alterations]. / Tabac et méthylation de l'ADN : un exemple d'altération épigénétique.

TOBACCO AND DNA METHYLATION: THE CASE FOR EPIGENETIC ALTERATIONS The mechanisms of the long-term impacts of exposure to chemical substances remain poorly understood. While genotoxic and mutagenic effects have been well characterized, epigenetic mechanisms such as DNA methylation could also account for the delayed effects of exposures. It is in the case of tobacco that the strongest arguments for a role of these mechanisms have been obtained in human populations. This text presents recent data on this issue demonstrating the plausibility of epigenetic mechanisms to explain the persistence of biological signals long after stopping exposure.

TABAC ET MÉTHYLATION DE L'ADN : UN EXEMPLE D'ALTÉRATION ÉPIGÉNÉTIQUE Les mécanismes des impacts à long terme des expositions à des substances chimiques demeurent assez mal connus. Si les effets génotoxiques et mutagènes ont été bien caractérisés, les mécanismes épigénétiques, comme la méthylation de l'ADN, pourraient aussi rendre compte des effets différés des expositions : dans les populations humaines, les arguments les plus solides ont été obtenus dans le cas du tabac. Les données récentes sur ce sujet démontrent la plausibilité des mécanismes épigénétiques pour expliquer la persistance de signaux biologiques longtemps après l'arrêt de l'exposition.

Associations between specific IgE sensitization to 26 respiratory allergen molecules and HLA class II alleles in the EGEA cohort.

BACKGROUND: Allergy, the most frequent immune disorder affecting 30% of the world's population, is the consequence of immunoglobin E (IgE) sensitization to allergens. Among the genetic factors suspected to be involved in allergy, the HLA class-II genomic region is a strong candidate. OBJECTIVE: To assess the association between HLA class-II alleles and specific IgE (sIgE) sensitization to a large number of respiratory allergen molecules. METHODS: The analysis relied on 927 participants of the EGEA cohort, including 497 asthmatics. The study focuses on 26 aeroallergens recognized by sIgE in at least 5% of the study population (determined with the MEDALL chip with sIgE ≥ 0.3 ISU) and 23 imputed HLA class-II alleles. For each sIgE sensitization and HLA class-II allele, we fitted a logistic regression model accounting for familial dependence and adjusted for gender, age, and genetic principal components. p-values were corrected for multiple comparisons (False Discovery Rate). RESULTS: Most of the 19 statistically significant associations observed regard pollen allergens (mugwort Art v 1, olive tree Ole e 1, timothy grass Phl p 2, Phl p 5 and plantain Pla l 1), three were mold allergen (Alternaria Alt a 1), and a single one regards house dust mite allergen (Der p 7). No association was observed with pet allergens. The strongest associations were found with mugwort Art v 1 (OR = 5.42 (95%CI, 3.30; 8.88), 4.14 (2.65; 6.47), 3.16 (1.88; 5.31) with DQB1*05:01, DQA1*01:01 and DRB1*01:01, respectively). CONCLUSION: Our results support the important role of HLA class-II alleles as immune response genes predisposing their carriers for sensitization to various major pollen allergens.

Keratoconus Prevalence in Families: A French Study.

PURPOSE: To assess the keratoconus (KC) prevalence among first-degree relatives of affected subjects. METHODS: We performed large epidemiological, prospective, and observational cohort, which took place in the French KC National Reference Centre, in Bordeaux. We studied data from 94 unrelated patients with KC and at least 2 first-degree relatives assessable. Each subject had a thorough ocular examination including visual acuity, refraction, slitlamp biomicroscopy, indirect ophthalmoscopy, keratometry, corneal topography and tomographer (Topographic Modeling System; Tomey Corporation, Japan; and Galilei G4, Ziemer, Switzerland, respectively), pachymetry, and biomechanical properties of cornea (Ocular Response Analyzer; Reichert, Depew, NY). Additional data were collected about atopy, smoking, and eye-rubbing habits. KC prevalence, multivariate analysis, and familial correlation analysis have been studied. RESULTS: Of the 221 relatives, 9.05% (n = 20) had clinically manifest KC, and 15.4% (n = 31) had subclinical KC. The prevalence of KC was estimated equal to 0.14 [95% confidence interval (CI), 0.07-0.22] among parents, 0.03 (95% CI, 0-0.10) among offspring, and 0.10 (95% CI, 0.04-0.20) among siblings. In multivariable analysis, age and eye-rubbing were associated with an increased risk of KC. The familial correlation of KC among probands' first-degree relatives was estimated at 0.55 among parents, 0.29 among offspring, and 0.49 among siblings. CONCLUSIONS: Being a first-degree relative of a subject with KC remains a risk factor for developing KC. This should be taken into account especially when screening for KC in pediatric patients, for refractive surgery candidates, and for patients being evaluated for cataract surgery.

Genome-wide interaction study of early-life smoking exposure on time-to-asthma onset in childhood.

BACKGROUND: Asthma, a heterogeneous disease with variable age of onset, results from the interplay between genetic and environmental factors. Early-life tobacco smoke (ELTS) exposure is a major asthma risk factor. Only a few genetic loci have been reported to interact with ELTS exposure in asthma. OBJECTIVE: Our aim was to identify new loci interacting with ELTS exposure on time-to-asthma onset (TAO) in childhood. METHODS: We conducted genome-wide interaction analyses of ELTS exposure on time-to-asthma onset in childhood in five European-ancestry studies (totalling 8273 subjects) using Cox proportional-hazard model. The results of all five genome-wide analyses were meta-analysed. RESULTS: The 13q21 locus showed genome-wide significant interaction with ELTS exposure (P = 4.3 × 10-8 for rs7334050 within KLHL1 with consistent results across the five studies). Suggestive interactions (P < 5 × 10-6 ) were found at three other loci: 20p12 (rs13037508 within MACROD2; P = 4.9 × 10-7 ), 14q22 (rs7493885 near NIN; P = 2.9 × 10-6 ) and 2p22 (rs232542 near CYP1B1; P = 4.1 × 10-6 ). Functional annotations and the literature showed that the lead SNPs at these four loci influence DNA methylation in the blood and are located nearby CpG sites reported to be associated with exposure to tobacco smoke components, which strongly support our findings. CONCLUSIONS AND CLINICAL RELEVANCE: We identified novel candidate genes interacting with ELTS exposure on time-to-asthma onset in childhood. These genes have plausible biological relevance related to tobacco smoke exposure. Further epigenetic and functional studies are needed to confirm these findings and to shed light on the underlying mechanisms.

Correction to: Functional variation in allelic methylomes underscores a strong genetic contribution and reveals novel epigenetic alterations in the human epigenome.

Following publication of the original article [1], the authors reported an error in Additional file 1.

Interactive effect between ATPase-related genes and early-life tobacco smoke exposure on bronchial hyper-responsiveness detected in asthma-ascertained families.

BACKGROUND: A positional cloning study of bronchial hyper-responsiveness (BHR) at the 17p11 locus in the French Epidemiological study on the Genetics and Environment of Asthma (EGEA) families showed significant interaction between early-life environmental tobacco smoke (ETS) exposure and genetic variants located in DNAH9. This gene encodes the heavy chain subunit of axonemal dynein, which is involved with ATP in the motile cilia function.Our goal was to identify genetic variants at other genes interacting with ETS in BHR by investigating all genes belonging to the 'ATP-binding' and 'ATPase activity' pathways which include DNAH9, are targets of cigarette smoke and play a crucial role in the airway inflammation. METHODS: Family-based interaction tests between ETS-exposed and unexposed BHR siblings were conducted in 388 EGEA families. Twenty single-nucleotide polymorphisms (SNP) showing interaction signals (p≤5.10-3) were tested in the 253 Saguenay-Lac-Saint-Jean (SLSJ) families. RESULTS: One of these SNPs was significantly replicated for interaction with ETS in SLSJ families (p=0.003). Another SNP reached the significance threshold after correction for multiple testing in the combined analysis of the two samples (p=10-5). Results were confirmed using both a robust log-linear test and a gene-based interaction test. CONCLUSION: The SNPs showing interaction with ETS belong to the ATP8A1 and ABCA1 genes, which play a role in the maintenance of asymmetry and homeostasis of lung membrane lipids.

A novel role for ciliary function in atopy: ADGRV1 and DNAH5 interactions.

BACKGROUND: Atopy, an endotype underlying allergic diseases, has a substantial genetic component. OBJECTIVE: Our goal was to identify novel genes associated with atopy in asthma-ascertained families. METHODS: We implemented a 3-step analysis strategy in 3 data sets: the Epidemiological Study on the Genetics and Environment of Asthma (EGEA) data set (1660 subjects), the Saguenay-Lac-Saint-Jean study data set (1138 subjects), and the Medical Research Council (MRC) data set (446 subjects). This strategy included a single nucleotide polymorphism (SNP) genome-wide association study (GWAS), the selection of related gene pairs based on statistical filtering of GWAS results, and text-mining filtering using Gene Relationships Across Implicated Loci and SNP-SNP interaction analysis of selected gene pairs. RESULTS: We identified the 5q14 locus, harboring the adhesion G protein-coupled receptor V1 (ADGRV1) gene, which showed genome-wide significant association with atopy (rs4916831, meta-analysis P value = 6.8 × 10-9). Statistical filtering of GWAS results followed by text-mining filtering revealed relationships between ADGRV1 and 3 genes showing suggestive association with atopy (P ≤ 10-4). SNP-SNP interaction analysis between ADGRV1 and these 3 genes showed significant interaction between ADGRV1 rs17554723 and 2 correlated SNPs (rs2134256 and rs1354187) within the dynein axonemal heavy chain 5 (DNAH5) gene (Pmeta-int = 3.6 × 10-5 and 6.1 × 10-5, which met the multiple-testing corrected threshold of 7.3 × 10-5). Further conditional analysis indicated that rs2134256 alone accounted for the interaction signal with rs17554723. CONCLUSION: Because both DNAH5 and ADGRV1 contribute to ciliary function, this study suggests that ciliary dysfunction might represent a novel mechanism underlying atopy. Combining GWAS and epistasis analysis driven by statistical and knowledge-based evidence represents a promising approach for identifying new genes involved in complex traits.

Multiancestry association study identifies new asthma risk loci that colocalize with immune-cell enhancer marks.

We examined common variation in asthma risk by conducting a meta-analysis of worldwide asthma genome-wide association studies (23,948 asthma cases, 118,538 controls) of individuals from ethnically diverse populations. We identified five new asthma loci, found two new associations at two known asthma loci, established asthma associations at two loci previously implicated in the comorbidity of asthma plus hay fever, and confirmed nine known loci. Investigation of pleiotropy showed large overlaps in genetic variants with autoimmune and inflammatory diseases. The enrichment in enhancer marks at asthma risk loci, especially in immune cells, suggested a major role of these loci in the regulation of immunologically related mechanisms.

Vitamin D levels and susceptibility to asthma, elevated immunoglobulin E levels, and atopic dermatitis: A Mendelian randomization study.

BACKGROUND: Low circulating vitamin D levels have been associated with risk of asthma, atopic dermatitis, and elevated total immunoglobulin E (IgE). These epidemiological associations, if true, would have public health importance, since vitamin D insufficiency is common and correctable. METHODS AND FINDINGS: We aimed to test whether genetically lowered vitamin D levels were associated with risk of asthma, atopic dermatitis, or elevated serum IgE levels, using Mendelian randomization (MR) methodology to control bias owing to confounding and reverse causation. The study employed data from the UK Biobank resource and from the SUNLIGHT, GABRIEL and EAGLE eczema consortia. Using four single-nucleotide polymorphisms (SNPs) strongly associated with 25-hydroxyvitamin D (25OHD) levels in 33,996 individuals, we conducted MR studies to estimate the effect of lowered 25OHD on the risk of asthma (n = 146,761), childhood onset asthma (n = 15,008), atopic dermatitis (n = 40,835), and elevated IgE level (n = 12,853) and tested MR assumptions in sensitivity analyses. None of the four 25OHD-lowering alleles were associated with asthma, atopic dermatitis, or elevated IgE levels (p ≥ 0.2). The MR odds ratio per standard deviation decrease in log-transformed 25OHD was 1.03 (95% confidence interval [CI] 0.90-1.19, p = 0.63) for asthma, 0.95 (95% CI 0.69-1.31, p = 0.76) for childhood-onset asthma, and 1.12 (95% CI 0.92-1.37, p = 0.27) for atopic dermatitis, and the effect size on log-transformed IgE levels was -0.40 (95% CI -1.65 to 0.85, p = 0.54). These results persisted in sensitivity analyses assessing population stratification and pleiotropy and vitamin D synthesis and metabolism pathways. The main limitations of this study are that the findings do not exclude an association between the studied outcomes and 1,25-dihydoxyvitamin D, the active form of vitamin D, the study was underpowered to detect effects smaller than an OR of 1.33 for childhood asthma, and the analyses were restricted to white populations of European ancestry. This research has been conducted using the UK Biobank Resource and data from the SUNLIGHT, GABRIEL and EAGLE Eczema consortia. CONCLUSIONS: In this study, we found no evidence that genetically determined reduction in 25OHD levels conferred an increased risk of asthma, atopic dermatitis, or elevated total serum IgE, suggesting that efforts to increase vitamin D are unlikely to reduce risks of atopic disease.

Functional variation in allelic methylomes underscores a strong genetic contribution and reveals novel epigenetic alterations in the human epigenome.

BACKGROUND: The functional impact of genetic variation has been extensively surveyed, revealing that genetic changes correlated to phenotypes lie mostly in non-coding genomic regions. Studies have linked allele-specific genetic changes to gene expression, DNA methylation, and histone marks but these investigations have only been carried out in a limited set of samples. RESULTS: We describe a large-scale coordinated study of allelic and non-allelic effects on DNA methylation, histone mark deposition, and gene expression, detecting the interrelations between epigenetic and functional features at unprecedented resolution. We use information from whole genome and targeted bisulfite sequencing from 910 samples to perform genotype-dependent analyses of allele-specific methylation (ASM) and non-allelic methylation (mQTL). In addition, we introduce a novel genotype-independent test to detect methylation imbalance between chromosomes. Of the ~2.2 million CpGs tested for ASM, mQTL, and genotype-independent effects, we identify ~32% as being genetically regulated (ASM or mQTL) and ~14% as being putatively epigenetically regulated. We also show that epigenetically driven effects are strongly enriched in repressed regions and near transcription start sites, whereas the genetically regulated CpGs are enriched in enhancers. Known imprinted regions are enriched among epigenetically regulated loci, but we also observe several novel genomic regions (e.g., HOX genes) as being epigenetically regulated. Finally, we use our ASM datasets for functional interpretation of disease-associated loci and show the advantage of utilizing naïve T cells for understanding autoimmune diseases. CONCLUSIONS: Our rich catalogue of haploid methylomes across multiple tissues will allow validation of epigenome association studies and exploration of new biological models for allelic exclusion in the human genome.

SigMod: an exact and efficient method to identify a strongly interconnected disease-associated module in a gene network.

MOTIVATION: Apart from single marker-based tests classically used in genome-wide association studies (GWAS), network-assisted analysis has become a promising approach to identify a set of genes associated with disease. To date, most network-assisted methods aim at finding genes connected in a background network, whatever the density or strength of their connections. This can hamper the findings as sparse connections are non-robust against noise from either the GWAS results or the network resource. RESULTS: We present SigMod, a novel and efficient method integrating GWAS results and gene network to identify a strongly interconnected gene module enriched in high association signals. Our method is formulated as a binary quadratic optimization problem, which can be solved exactly through graph min-cut algorithms. Compared to existing methods, SigMod has several desirable properties: (i) edge weights quantifying confidence of connections between genes are taken into account, (ii) the selection path can be computed rapidly, (iii) the identified gene module is strongly interconnected, hence includes genes of high functional relevance, and (iv) the method is robust against noise from either the GWAS results or the network resource. We applied SigMod to both simulated and real data. It was found to outperform state-of-the-art network-assisted methods in identifying disease-associated genes. When SigMod was applied to childhood-onset asthma GWAS results, it successfully identified a gene module enriched in consistently high association signals and made of functionally related genes that are biologically relevant for asthma. AVAILABILITY AND IMPLEMENTATION: An R package SigMod is available at: https://github.com/YuanlongLiu/SigMod. CONTACT: yuanlong.liu@inserm.fr. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.